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Abstract 

This  thesis  uses  statistical  analysis  to  forecast  the  probability  of  meeting  or 
exceeding  the  maximum  allowable  wind  speeds  for  each  of  the  launch  pads  at  the 
Kennedy  Space  Center  (KSC)  and  Cape  Canaveral  Air  Station  (CCAS).  Wind  data 
were  collected  from  the  Weather  Information  Network  Display  System  (WINDS), 
a  collection  of  47  meteorological  towers  located  throughout  KSC  and  CCAS,  over 
a  period  of  five  winters.  A  Fortran  program  was  written  to  calculate  conditional 
probabilities  of  meeting  or  exceeding  a  given  threshold  speed  during  eight  consecutive 
one-hour  periods,  using  the  current  wind  direction  and  peak  wind  speed  as  inputs. 
Forecast  probabilities  were  displayed  in  a  table  according  to  time  period  and  wind 
direction. 

Accuracy  was  measured  by  constructing  contingency  tables  and  calculating 
various  measures  of  accuracy.  Results  were  tested  for  significance  by  calculating 
p- values  for  the  chi-square  test.  This  method  was  found  to  have  very  little  skill  in 
forecasting  maximum  wind  speeds.  It  is  not  recommended  for  operational  use. 
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Use  of  Climatology  to  Predict  Maximum  Wintertime  Wind  Speeds  at 
the  Kennedy  Space  Center  and  Cape  Canaveral  Air  Station 


I.  Introduction 

1.1  Background 

Launch  Weather  Officers  (LWOs)  at  the  45th  Weather  Squadron  (45  WS), 
Patrick  Air  Force  Base  (PAFB),  are  responsible  for  making  “go/no-go”  decisions  dur¬ 
ing  the  space  launch  countdown  based  on  the  probability  of  violating  various  weather 
constraints.  One  of  the  biggest  problems  the  LWOs  face  is  accurately  forecasting 
peak  winds  during  the  winter  season,  which  is  defined  as  November  through  March,  at 
the  Kennedy  Space  Center  (KSC)  and  Cape  Canaveral  Air  Station  (CCAS).  About 
eight  hours  before  the  launch  time,  the  structure  that  supports  the  rocket  and  pro¬ 
tects  it  from  the  elements  is  rolled  away  from  the  rocket,  leaving  it  exposed  and 
unsupported.  If  the  wind  speeds  are  too  strong  during  this  eight-hour  countdown,  a 
collision  could  occur  between  the  rocket  and  its  supporting  structure.  Therefore,  if 
the  wind  speeds  are  expected  to  exceed  a  particular  threshold,  which  varies  depend¬ 
ing  on  the  launch  vehicle  and  the  wind  direction,  the  launch  must  either  be  delayed 
or  rescheduled.  Of  the  32  NASA  and  Air  Force  launches  in  1999,  six  were  delayed 
due  to  weather-related  problems,  costing  the  government  up  to  $616,000  per  day  to 
postpone  a  launch  (Frank,  2000  and  White,  2000). 

Strong  winds  at  KSC  and  CCAS  can  occur  at  any  time  of  the  year.  In  the 
summer  season,  these  winds  are  usually  caused  by  thunderstorms  and  microbursts. 
Thunderstorms  are  rare  in  the  winter  season,  however,  and  frontal  passages  are  the 
main  cause  of  strong  winds.  The  four  types  of  fronts  affecting  central  Florida  are 
discussed  in  Chapter  2. 
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Wind  data  are  collected  using  one  of  the  most  densely  instrumented  meteoro¬ 
logical  mesonetworks  in  the  United  States,  the  Weather  Information  Network  Display 
System  (WINDS).  It  consists  of  47  meteorological  towers  located  throughout  KSC 
and  CCAS  and  covers  an  area  of  approximately  1200  km2,  with  an  average  station 
density  of  one  tower  per  27  km2.  All  of  these  towers  measure  wind  speed,  wind  di¬ 
rection,  and  temperature,  and  selected  towers  also  measure  relative  humidity.  The 
sensors  on  these  towers  range  in  height  from  6  feet  to  492  feet.  Observations  are 
recorded  at  one-  and  five-minute  intervals  for  a  total  of  up  to  8580  observations  each 
minute  (Raytheon,  1998). 

1.2  Problem  Statement 

Personnel  at  the  45  WS  have  a  vast  amount  of  information  available  to  them 
via  the  WINDS  network.  However,  the  LWOs  have  little  guidance  to  help  them 
forecast  maximum  wind  speeds  during  the  winter.  Storch  (1999)  attempted  to 
forecast  wintertime  wind  speeds  for  KSC  and  CCAS  using  a  neural  network,  but  his 
results  were  inconclusive.  In  addition,  there  have  been  no  previous  climatological 
studies  of  wintertime  wind  speeds  for  KSC  and  CCAS.  Thus,  the  goal  of  this  thesis 
is  to  design  a  climatology-based,  short-term  (less  than  8  horns)  forecasting  tool  to 
assist  the  LWOs  with  their  wind  forecasts. 

1.3  Scope 

This  project  focuses  on  developing  a  conditional  climatology  program  to  help 
the  LWOs  predict  peak  winds  more  accurately.  Personnel  at  45  WS  requested  that 
only  the  WINDS  towers  located  close  to  the  launch  pads  be  used  in  this  study.  These 
towers  are  listed  in  Table  1  and  can  be  found  on  the  map  presented  in  Figure  1. 

Only  data  from  the  1994-95,  1995-96,  1996-97,  1997-98,  and  1998-99  winter 
seasons  have  been  included  in  this  study  because  unreliable  equipment  reported 
inaccurate  wind  observations  at  the  launch  pads  prior  to  1994  (Roeder,  1999).  In 
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Tower 

Level  (feet) 

Side 

Location 

Rocket 

NW 

Pad  17 

SE 

Pad  17 

0036 

90 

M 

lESElHl 

Atlas 

mm 

54 

Titan 

noi 

162 

NW 

Pad  40 

Titan 

mm 

54 

■SMI 

Titan 

fBI 

162 

Titan 

E 

NW 

Shuttle 

SE 

Shuttle 

0393 

60 

NW 

Pad  39B 

Shuttle 

0394 

60 

SE 

Pad  39B 

Shuttle 

Table  1  WINDS  towers  used  in  this  study. 


addition,  previous  studies  on  mesoscale  systems  such  as  thunderstorms,  microbursts, 
and  sea  breezes  have  vastly  improved  wind  forecasts  during  the  summer  (Roeder, 
1999).  The  winter  season  continues  to  be  a  problem  because  LWOs  are  unfamiliar 
with  how  the  relatively  infrequent  frontal  passages  affect  wind  speed,  especially  when 
the  peak  winds  are  already  within  a  few  knots  of  the  threshold  value. 


1-4  Thesis  Organization 

The  remainder  of  this  thesis  explores  the  viability  of  predicting  maximum  wind 
speeds  over  an  eight-hour  period  at  KSC  and  adjacent  CCAS  using  conditional  cli¬ 
matology  developed  from  a  limited  database.  Chapter  2  presents  a  brief  overview  of 
weather  systems  which  affect  central  Florida  and  discussed  recent  studies  which  have 
looked  at  this  same  problem.  Chapter  3  includes  sensor  accuracy  and  a  description 
of  the  data,  including  a  discussion  of  observed  winds  and  their  causes,  as  well  as 
statistical  methods  used  in  the  analysis  of  this  data.  Chapter  4  describes  the  steps 
taken  to  conduct  this  study  and  discusses  the  results.  Chapter  5  draws  conclusions 
about  how  the  algorithm  performed.  Observed  wind  frequency  tables  for  all  of  the 
WINDS  towers  used  in  this  study  and  contingency  tables  used  in  the  verification  of 
the  output  are  contained  in  the  appendices. 
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Figure  1  Location  of  the  WINDS  towers.  Dark  gray  triangles  represent  WINDS 
towers.  Thin  rectangles  represent  launch  pads.  Towers  of  interest  are 
0002,  0036,  0393,  0394,  0397,  0398,  1101,  and  1102. 
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II.  Literature  Review 


Due  to  their  locations  in  the  northern  hemisphere  subtropics,  PAFB  and  CCAS 
are  affected  by  the  subtropical  ridge  known  as  the  Bermuda  high  throughout  the 
year.  The  central  Florida  coast  experiences  two  main  seasons:  the  airmass  regime, 
which  runs  from  mid-May  to  late  September,  and  the  frontal  regime,  which  takes 
place  from  late  November  to  mid-March.  During  the  frontal  regime,  winds  can  be 
tricky  to  forecast  as  relatively  infrequent  frontal  passages  occur.  Additionally,  two 
short  transition  seasons  take  place,  one  in  the  spring  and  another  in  the  fall.  A 
brief  overview  of  each  of  the  seasons  is  presented  here.  Seasonal  descriptions  were 
obtained  from  the  Terminal  Forecast  Reference  Notebook  (1987). 

2.1  Airmass  Regime 

The  Bermuda  high  is  strongest  during  the  airmass  regime,  when  it  is  centered 
around  30  degrees  north  latitude.  The  position  and  strength  of  the  high  prevents 
frontal  passages  from  advancing  over  Florida,  as  shown  in  Figure  2.  This  season  is 
dominated  by  maritime  tropical  air,  also  due  to  the  position  of  the  high.  The  days 
are  characterized  by  a  nearly  continuous  easterly  wind  and  convective  activity,  both 
of  which  can  cause  strong  winds. 

2.2  Frontal  Regime 

Florida’s  frontal  regime  takes  place  from  late  November  to  mid-March.  This 
regime  is  a  mild  winter,  during  which  the  Bermuda  high  weakens  and  retreats  to 
around  20°  north  latitude,  and  true  frontal  passages  occur.  The  four  types  of 
cold  fronts  which  affect  central  Florida  during  the  winter  are  categorized  accord¬ 
ing  to  their  orientations:  north/south,  northeast /southwest,  east/west,  and  north¬ 
west/southeast.  All  of  these  fronts  can  cause  strong  winds  at  the  launch  pads. 


5 


Figure  2  Position  of  the  Bermuda  high  during  the  summer  season. 

2.2.1  North/South  Oriented  Fronts.  The  north/south  oriented  front  is 
associated  with  a  maritime  polar  front.  As  shown  in  Figure  3,  the  low-pressure 
center  is  positioned  over  east-central  Alabama  or  west-central  Georgia,  and  the  high- 
pressure  center  is  located  over  or  just  north  of  the  central  Gulf  of  Mexico.  This 
synoptic  situation  can  cause  surface  winds  in  excess  of  35  knots  along  the  east  coast 
of  central  Florida. 

2.2.2  Northeast/Southwest  Oriented  Fronts.  The  northeast/southwest  ori¬ 
ented  front  can  occur  with  either  a  maritime  polar  or  continental  polar  airmass.  The 
low-pressure  center  is  usually  located  between  North  Carolina  and  Bermuda  and  the 
high-pressure  center  is  located  over  Arkansas  or  Tennessee.  “This  type  of  front 
produces  the  most  severe  weather  in  Florida  in  the  form  of  thunderstorms  with  high 
winds  and  hail”  ( Terminal  Forecast  Reference  Notebook,  1987).  However,  because 
the  front  usually  decelerates  and  weakens  when  it  reaches  the  Florida  Panhandle, 
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most  of  the  severe  weather  occurs  in  the  northern  part  of  the  state.  See  Figure  4 
for  an  example  of  this  situation. 

2.2.3  East/ West  Oriented  Fronts.  East/west  oriented  fronts  are  associated 
with  a  continental  polar  airmass.  The  low-pressure  center  is  located  northeast  of 
Bermuda  and  the  high-pressure  center  is  located  in  eastern  Ohio  or  western  Penn¬ 
sylvania.  Northeasterly  winds  behind  the  front  can  exceed  20  knots.  See  Figure  5 
for  an  example  of  this  situation. 

2.2.4  Northwest/Southeast  Oriented  Fronts.  An  east/west  oriented  front 
occasionally  stalls  between  Jacksonville  and  Daytona  Beach.  Sometimes  it  begins 
to  move  south  again,  this  time  with  a  northwest/southeast  orientation  as  shown  in 
Figure  6.  Strong  westerly  winds  occur  ahead  of  the  front  and  become  northeasterly 
behind  the  front. 

2.3  Transition  Seasons 

The  spring  transition  season  lasts  from  late  March  to  mid-May,  when  the  sub¬ 
tropical  ridge  strengthens  and  shifts  to  the  north.  Frontal  passages  are  less  frequent, 
and  winds  become  easterly. 

The  fall  transition  season  lasts  from  early  October  to  mid-November.  Weak 
fronts  begin  to  move  through  the  area,  and  winds  become  northerly. 

2/  Previous  Studies 

As  mentioned  in  Chapter  1,  the  problem  of  predicting  winds  at  KSC/CCAS 
was  examined  by  Storch  (1999),  who  developed  a  neural  network  to  forecast  winds 
for  the  launch  pads.  Storch’s  results  showed  that  his  neural  network  performed 
no  better  than  persistence.  Cloys  (2000)  continued  research  in  the  neural  network 
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field.  Cloys  developed  his  network  using  a  different  method  and  found  that  his 
model  performed  better  than  persistence  late  in  the  forecast  period. 
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Bermuda 


Figure  4  Typical  synoptic  pattern  for  the  Northeast/Southwest  oriented  front 


Figure  6  Typical  synoptic  pattern  for  the  Northwest /Southeast  oriented  front. 


III.  Data  Description  and  Analysis 


3.1  Data  Description 

Personnel  at  the  45  WS  provided  wind  observations  gathered  via  the  WINDS 
network  from  the  1994-95,  1995-96,  1996-97,  1997-98,  and  1998-99  winter  seasons. 
Observations  were  taken  at  five-minute  intervals.  Over  the  course  of  five  winters, 
this  adds  up  to  217728  observations  for  each  tower  and  level,  assuming  no  data  are 
missing  data. 

3.1.1  Sensor  Accuracy.  The  WINDS  towers  used  in  this  study  are  equipped 
with  the  R.  M.  Young  Model  05305-18  Wind  Monitor- AQ  sensor  (Raytheon,  1998). 
These  sensors  are  accurate  to  within  0.3  ms-1  for  speed  and  within  3°  for  direction 
Other  sensor  characteristics  are  listed  in  Table  2. 

3.1.2  Description  of  Observed  Data.  Tables  3  through  8  contain  observed 
wind  speed  and  direction  frequencies  for  Tower  1101  at  the  162-foot  level.  Al¬ 
though  wintertime  winds  can  occur  from  any  direction,  two  sectors  appear  to  be 
predominant  throughout  the  winter  although  there  is  a  slight  variation  from  month 
to  month.  As  shown  in  Table  3,  the  most  frequently-appearing  wind  directions  are 
west  through  north,  which  occur  in  conjunction  with  frontal  passages.  A  smaller  but 
still  significant  maximum  is  seen  from  the  southeast  through  south,  which  is  some¬ 
times  accounted  for  by  the  Bermuda  high  and  sometimes  caused  by  the  approach  of 
a  front. 


Characteristic 

Speed  (ms  x) 

Direction  (degrees) 

Operating  range 

0.33  -  54 

0-360 

Survival  range 

60 

N/A 

Accuracy 

±0.3 

d=3 

Starting  Threshold 

0.33 

<0.33  for  10°  displacement 

Table  2  Characteristics  of  wind  speed  and  direction  sensors  (Raytheon.  1998). 
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In  November  (Table  4),  winds  are  most  frequently  from  300°  to  329°,  due 
to  frontal  passages,  and  0°  to  089°,  partly  due  to  frontal  passages  and  partly  due 
to  the  onshore  gradient  caused  by  the  Bermuda  high.  Winds  from  the  southeast 
through  west  are  relatively  infrequent.  By  December  (Table  5),  westerly  through 
northerly  winds  (270°  to  359°)  occur  with  a  much  greater  frequency  and  a  slight 
increase  in  southeasterly  and  southwesterly  winds  (150°  to  209°)  begins  to  appear. 
This  is  also  evident  in  January  (Table  6)  and  February  (Table  7).  However,  in 
February,  the  southerly  winds  become  more  southeasterly,  shifting  to  between  120° 
and  179°.  By  March  (Table  8),  frontal  passages  begin  to  occur  less  frequently  with 
the  strengthening  of  the  Bermuda  high.  Therefore,  northerly  winds  are  seen  less 
frequently  while  southeasterly  winds  become  more  common. 

Another  feature  worth  noting  is  the  frequency  of  wind  speed.  The  vast  ma¬ 
jority  of  winds  are  between  5  and  19  knots,  regardless  of  the  month.  Wind  speeds 
of  4  knots  or  less  generally  occur  less  than  7%  of  the  time.  On  average,  wind  speeds 
greater  than  or  equal  to  20  knots  occur  roughly  19%  of  the  time  and  winds  greater 
than  or  equal  to  25  knots  occur  about  7%  of  the  time.  Wind  speeds  greater  than 
40  knots  are  extremely  rare,  although  a  slight  increase  in  the  frequency  of  stronger 
winds  occurs  from  January  to  February. 

Observed  frequency  tables  for  the  other  WINDS  towers  used  in  this  study  can 
be  found  in  Appendix  A. 

3.2  Data  Analysis 

3.2.1  Conditional  Climatology.  Conditional  climatology  is  frequently  used 
by  the  meteorological  community  to  quantify  the  likelihood  that  a  particular  event 
will  occur,  such  as  visibility  decreasing  to  below  one-half  mile  within  the  next  hour, 
given  that  the  visibility  is  currently  two  miles.  Two  specific  examples  of  conditional 
climatology  products  commonly  used  in  an  Air  Force  weather  station  are  Modeled 
Ceiling  and  Visibility  (MODCV)  and  Wind-Stratified  Conditional  Climatology  Ta- 
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10-14 

15-19 

20-24 

25-29 

30-34 

35-39 

>40 

Total 

0-29 

0.33 

1.97 

wmm 

2.57 

1.24 

0.06 

0.02 

iilmi 

30-59 

0.22 

1.77 

1.20 

1.15 

0.60 

0.10 

0.11 

60-89 

0.26 

1.87 

1.78 

2.23 

1.41 

0.41 

0.10 

0.06 

0.01 

90-119 

0.28 

2.51 

2.44 

2.03 

0 

mm 

0 

8.15 

120-149 

0.33 

2.47 

2.84 

2.93 

WBfiM 

0.48 

0.04 

0 

10.6 

150-179 

2.69 

HIM 

0.61 

0.28 

9.25 

180-209 

1.85 

1.92 

0.21 

0.02 

7.30 

210-239 

0.27 

1.87 

1.58 

1.35 

0.63 

0.24 

0.04 

6.00 

240-259 

0.29 

1.97 

EIOM 

lisa 

0.62 

0.36 

0.13 

E2B 

6.88 

270-299 

0.29 

1.94 

2.71 

wmm 

1.19 

0.85 

0.24 

mm 

ISIiEl 

9.53 

0.31 

2.59 

2.90 

2.13 

1.32 

mm 

0 

9.55 

0.22 

2.21 

2.57 

2.32 

1.82 

mm 

0.07 

10.4 

Total 

26.7 

24.5 

0.55 

0.19 

100 

Table  3  Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
1101,  Level  162,  based  on  130314  observations.  Wind  speeds  (first  row) 
are  in  knots  and  wind  directions  (first  column)  are  in  degrees. 


10-14 

15-19 

35-39 

>40 

Total 

0-29 

2.22 

4.00 

0 

0 

11.1 

30-59 

0.30 

1.58 

1.40 

2.10 

2.42 

1.11 

0.42 

0.53 

0.03 

9.9 

EEH 

3.18 

0.51 

McM 

mm 

pm 

1.48 

mm 

IBlB 

0.40 

2.62 

0.40 

0 

0 

0 

0.39 

2.10 

0.77 

PW1 

0 

0 

0 

180-209 

0.10 

1.92 

1.18 

0.05 

0 

0 

0 

210-239 

WEEl 

mm 

Eml 

0.07 

0 

0 

Kin 

240-259 

EMI 

mm 

mm 

0.03 

mm 

0 

0 

HM 

0.59 

2.02 

m<w 

0.38 

0 

0 

KfH^l 

0.69 

2.80 

1.16 

BKT31 

0 

0 

330-359 

0.33 

1.71 

2.59 

1.68 

0.57 

0.38 

0 

0 

0 

7.27 

Total 

4.94 

26.1 

24.9 

24.5 

13.6 

4.06 

0.94 

0.87 

0.03 

100 

Table  4  Percent  observed  wind  speed  and  directions  for  November,  Tower  1101, 
Level  162,  based  on  17695  observations.  Wind  speeds  (first  row)  are  in 
knots  and  wind  directions  (first  column)  are  in  degrees. 
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Table  5  Percent  observed  wind  speed  and  directions  for  December,  Tower  1101, 
Level  162,  based  on  23888  observations.  Wind  speeds  (first  row)  are  in 
knots  and  wind  directions  (first  column)  are  in  degrees. 


0-4 

5-9 

30-34 

35-39 

>40 

Total 

0-29 

0.26 

1.58 

0.02 

0 

0 

6.73 

30-59 

0.25 

2.15 

0.85 

1.55 

0.33 

0 

0 

0 

0 

5.13 

60-89 

0.25 

1.53 

1.68 

2.95 

1.31 

0 

7.76 

90-119 

0.39 

1.72 

2.12 

3.09 

0.32 

0 

7.64 

120-149 

0.54 

2.45 

I1M 

0.77 

0.05 

0 

0 

7.63 

150-179 

0.45 

1.99 

1.57 

0.62 

0.02 

0 

9.08 

180-209 

2.43 

1.41 

0.25 

0.01 

0 

0 

9.84 

210-239 

0.41 

0.33 

0.04 

0.01 

0 

7.05 

240-259 

0.33 

0.31 

0.17 

0.04 

0 

6.54 

270-299 

0.41 

2.24 

0.01 

11.9 

300-329 

0.26 

2.81 

3.82 

1 

mm 

0 

11.5 

330-359 

2.49 

1.79 

1.47 

0.69 

0.08 

9.08 

Total 

4.11 

25.9 

25.1 

25.9 

12.9 

4.78 

1.07 

0.15 

0.01 

100 

Table  6  Percent  observed  wind  speed  and  directions  for  January,  Tower  1101,  Level 
162,  based  on  22411  observations.  Wind  speeds  (first  row)  are  in  knots 
and  wind  directions  (first  column)  are  in  degrees. 
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0-29 
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1.70 

2.12 
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0.93 

0.21 

0 

0 

0 
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0.10 

1.42 

1.12 

0.93 

0.06 

0 

0 

0 

0 

3.63 
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0.12 

1.19 

0.28 

0.32 

0.06 

0 

0 

0 

0 

1.97 

90-119 

0.15 

1.96 

1.82 

1.29 

1.02 

0.44 

0 

0 

0 

6.68 

120-149 

0.15 

3.09 

2.68 

2.03 

1.49 

1.07 

0.16 

0.02 

0 

10.7 
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0.07 

2.33 

2.84 

2.92 

1.53 
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0.33 

0.30 

0.07 

11.5 
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2.54 

1.98 
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0.37 

0.06 

0 

0 

7.56 
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0.12 
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2.08 

1.86 

0.68 

0.06 

0.01 

0 

0 

6.81 

240-259 

0.17 

2.35 

2.42 

1.93 

0.62 

0.24 

0.16 

0.26 

0.11 

8.24 

270-299 

0.23 

2.56 

3.15 

1.91 

1.16 

1.21 

0.61 

0.60 

0.13 

11.6 

300-329 

0.35 

3.42 

2.80 

2.32 

1.70 

0.68 

0.11 

0.04 

0 

11.4 
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0.20 

2.67 

3.74 

3.59 

2.37 
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0.01 

0 

0 

13.3 

Total 

2.02 

26.5 

27.6 

22.6 

12.3 

6.02 

1.45 

1.22 

0.32 

100 

Table  7  Percent  observed  wind  speed  and  directions  for  February,  Tower  1101, 
Level  162,  based  on  22593  observations.  Wind  speeds  (first  row)  are  in 
knots  and  wind  directions  (first  column)  are  in  degrees. 
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Total 

0-29 

0.29 

2.14 

2.53 

2.80 

1.27 

0.71 

0.15 

0.06 

0.03 

9.98 

30-59 

0.24 

1.82 

1.33 

0.66 

0.27 

0.22 

0.07 

0.06 

0.03 

4.70 

60-89 

0.25 

2.21 

2.44 

2.46 

1.40 

0.25 

0.03 

0.02 

0.01 

9.05 

90-119 

0.24 

3.01 

3.29 

2.08 

0.67 

0.16 

0.01 

0 

0.01 

9.46 

120-149 

0.28 

2.09 

3.59 

4.54 

2.31 

0.59 

0.03 

0 

0 

13.4 

150-179 

0.16 

1.36 

3.12 

2.06 

1.40 

0.57 

0.34 

0.08 

0.01 
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0.20 

1.66 

1.62 

1.58 
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0.18 

0.02 

0 

0 

6.11 

210-239 

0.31 

1.75 

1.34 

1.00 

0.55 

0.35 

0.07 

0.03 

0 

5.39 

240-259 

0.28 

1.39 

2.22 

1.43 

1.01 

0.58 

0.14 

0.04 

0 

7.10 

270-299 

0.13 

1.43 

2.18 

2.54 

1.21 

0.57 

0.08 

0.02 

0 

8.17 

300-329 

0.15 

1.96 

2.15 

1.42 

0.96 

0.10 

0 

0 

0 

6.76 

330-359 

0.21 

2.08 

1.99 

2.20 

2.22 

1.47 

0.34 

0.08 

0 

10.7 

Total 

2.75 

22.9 

27.8 

24.8 

14.1 

5.76 

1.26 

0.38 

0.18 

100 

Table  8  Percent  observed  wind  speed  and  directions  for  March,  Tower  1101,  Level 
162,  based  on  43727  observations.  Wind  speeds  (first  row)  are  in  knots 
and  wind  directions  (first  column)  are  in  degrees. 
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SCT  (B) 

BKN  (B/) 

Horizontal  Totals 

TSRA  (A) 

35 

65 

100 

No  TSRA  (A/) 

65 

35 

Vertical  Totals 

100 

100 

Table  9  Conditional  probabilities  for  the  occurrence  of  thunderstorms  versus  cloud 
cover. 

bles  (WSCC),  both  of  which  are  produced  by  the  Air  Force  Combat  Climatology 
Center. 

As  early  as  1905,  meteorologists  were  producing  probabilistic  forecasts  based 
on  simple  contingency  tables  (Murphy  and  Winkler,  1984).  A  greater  interest  in 
evaluating  these  probabilistic  forecasts  did  not  come  about  until  the  mid-1950s,  with 
the  advances  in  computer  technology  and  an  increased  interest  in  applying  statistical 
techniques  to  produce  weather  forecasts  (Murphy  and  Winkler,  1984).  Some  of  these 
same  methods  are  still  used  today,  nearly  100  years  later,  and  are  described  in  the 
remainder  of  this  chapter. 

3.2.2  Conditional  Probability.  According  to  Wilks  (1995),  conditional 
probability  is  the  probability  that  a  specific  event  will  occur,  given  that  some  other 
event  has  occurred  or  will  occur  in  the  future.  The  “conditioning  event”  is  the  event 
upon  which  subsequent  events  are  contingent. 

In  general,  conditional  probability  can  be  expressed  by  P(A  \  B)  =  , 

where  B  is  the  conditioning  event  and  A  is  the  subsequent  event. 

Conditional  probability  can  be  explained  by  using  the  example  of  a  thunder¬ 
storm  (TSRA)  forecast  based  on  cloud  cover.  Suppose  an  amateur  meteorologist 
counts  the  number  of  partly  cloudy  (SCT)  days  with  and  without  thunderstorms, 
and  the  number  of  mostly  cloudy  (BKN)  days  with  and  without  thunderstorms.  The 
totals  are  then  entered  into  a  table  similar  to  Table  9. 

To  calculate  the  probability  of  thunderstorm  occurrence  given  that  it  is  now 
mostly  cloudy,  simply  divide  the  number  of  mostly  cloudy  days  with  thunderstorms, 
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( A  n  B'),  by  the  total  number  of  mostly  cloudy  days,  (£?').  In  this  case,  ^  =  .65 
or  65%.  Similarly,  the  probability  that  it  is  mostly  cloudy  given  that  there  are 
no  thunderstorms  can  be  computed  by  dividing  the  number  of  mostly  cloudy  days 
without  thunderstorms,  (B'CiA'),  by  the  total  number  of  days  without  thunderstorms 
(A/).  This  comes  out  to  ^  =  .35  or  35%. 

3.2.3  Contingency  Tables.  Contingency  tables  can  be  used  to  extract  mea¬ 
sures  of  accuracy  from  yes/no  forecasts.  According  to  Wilks  (1995),  a  contingency 
table  is  used  to  display  the  possible  combinations  of  forecasted  and  observed  event 
pairs.  An  example  of  a  2  x  2  contingency  table  is  shown  in  Figure  7.  A  perfect 
forecast  would  produce  zeroes  in  blocks  b  and  c,  the  total  number  of  “yes”  forecasts 
in  block  a,  and  the  total  number  of  “no”  forecasts  in  block  d.  Imperfect  forecasts 
need  to  be  measured  in  some  way.  Several  of  these  measures  are  summarized  in  the 
following  section. 
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c 

d 

dd 

ale 

m 

n=aft>icid 

Figure  7  Contingency  Table  According  to  Wilks  (1995). 
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3.2.4  Measures  of  Accuracy.  A  measure  of  accuracy  describes  the  average 
correspondence  between  a  forecast  event  and  the  event  that  actually  occurred  (Wilks, 
1995). 


3.2.4-  1  Hit  Rate  (H).  The  hit  rate  is  a  basic  measure  of  accuracy 
and  can  be  described  as  the  proportion  of  correct  forecasts.  Using  the  terminology 
introduced  in  Figure  6,  the  hit  rate  is  calculated  by: 

H  =  where  a  is  the  number  of  times  strong  winds  were  forecast  and 
observed,  d  is  the  number  of  times  strong  winds  were  not  forecast  and  not  observed, 
and  n  is  the  total  sample  size. 

3.2. 4- 2  Critical  Success  Index  (CSI).  Also  known  as  the  threat  score, 
the  critical  success  index  is  the  number  of  correct  “yes”  forecasts  divided  by  the 
number  of  times  strong  winds  were  forecast  and  observed.  The  critical  success 
index  is  given  by: 

CSI  =  g+°+c  where  b  is  the  number  of  times  strong  winds  were  forecast  but 
did  not  occur,  and  c  is  the  number  of  times  strong  winds  occurred  but  were  not 
forecast. 


3.2.4-S  Probability  of  Detection  (POD).  The  probability  of  detection 
is  the  likelihood  that  strong  winds  were  forecast,  given  that  strong  winds  occurred. 
The  probability  of  detection  is  given  by: 

POD  — 


3.2. 4-4  False  Alarm  Rate  (FAR).  The  false  alarm  rate  is  the  pro¬ 
portion  of  forecast  events  that  do  not  occur  (Wilks,  1995).  The  false  alarm  rate  is 
given  by: 


far  =  Jj. 
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3. 2-4- 5  Bias  (B).  Bias  is  the  comparison  of  the  average  forecast  with 
the  average  observation,  and  is  usually  expressed  as  a  ratio  (Wilks,  1995).  A  bias 
equal  to  1  indicates  that  strong  winds  were  forecast  the  same  number  of  times  they 
occurred.  A  bias  less  than  1  indicates  that  strong  winds  were  underforecast,  while  a 
bias  greater  than  1  indicates  that  strong  winds  were  overforecast.  Bias  is  computed 
by: 

r>  __  a+b 

CL-\~C 

3.2.5  Skill  Score.  A  skill  scores  is  a  measure  of  relative  accuracy.  One 
should  not  rely  on  skill  score  alone  because  different  skill  scores  perform  differently. 
Wilks  (1995)  explains  this  inconsistency  by  stating  that  skill  scores  are  “scalar  mea¬ 
sures  of  forecast  performance  in  what  is  intrinsically  a  higher-dimensional  setting.” 
The  most  popular  measure  of  forecast  accuracy  is  the  Heidke  skill  score. 

3.2.5. 1  Heidke  Skill  Score  (HSS).  The  Heidke  skill  score  uses  the  hit 
rate  as  the  basic  measure  of  accuracy.  The  reference  measure  of  accuracy  is  the 
hit  rate  that  would  be  achieved  by  random  forecasts,  provided  that  the  marginal 
distributions  of  forecasts  and  observations  in  the  contingency  table  for  the  random 
forecasts  is  equal  to  the  marginal  distributions  in  the  verification  data  set  (Wilks, 
1995).  Using  Figure  6  as  an  example,  ^  and  ~  in  the  forecast  set  should  be  the 
same  as  ^  and  —  in  the  verification  set.  A  perfect  forecast  has  a  Heidke  score 
equal  to  one.  The  Heidke  skill  score  is  computed  as  follows: 

tt  c  Q  _ 2(ad—bc) _ 

(a+c)(c+d)+(a+fr)(6+d) 

3.2.6  Testing  for  Significance.  Statistical  measures  of  accuracy,  bias,  and 
skill  are  meaningless  unless  a  connection  between  the  observations  (columns)  and 
forecasts  (rows)  in  the  contingency  table  can  be  shown.  Otherwise,  the  observed 
and  forecast  data  are  said  to  be  independent  and  seemingly-impressive  accuracy 
measures  are  obtained  by  chance.  Results  were  tested  for  significance  at  the  0.05 
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level  by  calculating  p-values  for  the  chi-square  test.  If  the  p-value  is  less  than  0.05, 
the  notion  of  independence  is  rejected  and  the  type  of  wind  speed  observed  (yes/no) 
is  considered  dependent  on  the  forecast  category. 
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IV.  Methodology 


4-1  Preparing  the  Observations 

Observations  from  the  WINDS  network  were  first  filtered  with  a  simple  Fortran 
program  which  removed  all  but  the  towers  and  levels  of  interest  listed  in  Table  1. 
Because  this  study  only  used  wind  direction  and  peak  speed  as  inputs  and  in  its 
calculations,  unnecessary  information  such  as  temperature  and  relative  humidity  was 
also  removed.  The  remaining  data  were  organized  into  separate  text  files  according 
to  month,  tower  number,  and  level.  For  example,  wind  data  for  tower  1101  at  the 
54-foot  level  from  February  1995,  February  1996,  February  1997,  and  February  1998 
were  collected  in  a  file  named  “febll01-54.txt.”  Another  file  called  “feb99-1101- 
54.txt”  was  used  to  store  data  for  tower  1101  at  the  54-foot  level  from  February 
1999.  Data  from  the  second  file  were  used  to  verify  the  conditional  probability 
forecasts  based  on  data  from  the  first  file.  After  the  forecasting  algorithm  was 
tested,  both  of  these  files  were  combined  for  operational  use,  assuming  output  from 
the  algorithm  would  prove  to  be  valuable. 

Erroneous  observations  were  removed  manually.  Plots  of  peak  speed  versus 
time  and  wind  direction  versus  time  were  created  for  each  month  (see  examples  in 
Figures  8  and  9).  Any  obvious  spikes  on  the  speed  plot  were  compared  to  speeds 
leading  up  to  and  following  the  spike,  as  well  as  direction  for  the  same  time.  If  the 
spike  was  preceded  by  or  led  up  to  strong  winds  and  was  accompanied  by  a  wind 
shift,  it  was  considered  to  be  correct  and  was  not  removed  from  the  file.  If  the 
spike  was  not  accompanied  by  strong  winds  or  a  wind  shift,  it  was  considered  to  be 
erroneous  and  removed  from  the  file.  Only  two  observations  were  removed  using 
this  method. 

Next,  the  data  were  screened  for  missing  values.  If  either  the  wind  direction, 
speed,  or  peak  speed  were  missing,  all  three  values  were  replaced  with  “-9999.”  This 
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Wind  direction  (degrees) 


Figure  9  Sample  of  wind  directions  from  Tower  1101,  162-foot  level,  February 
1999. 
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replacement  value  allowed  the  time  position  to  be  retained  in  the  dataset  while  telling 
the  Fortran  program  (discussed  in  the  following  section)  to  disregard  these  values  in 
its  calculations. 

4-2  Calculating  the  Conditional  Probabilities 

A  Fortran  program  was  written  to  calculate  the  conditional  probabilities  of 
meeting  or  exceeding  a  certain  speed,  hereafter  referred  to  as  the  threshold  speed, 
using  the  current  wind  direction,  current  peak  speed,  and  threshold  speed  (which 
varies  depending  on  wind  direction  and  type  of  rocket)  as  inputs.  These  forecast 
probabilities  were  created  for  eight  one-hour  periods  following  the  initial  time. 

First,  the  user  selects  the  current  month,  tower  number,  and  sensor  height. 
This  allows  the  program  to  open  the  appropriate  file.  Current  direction,  peak 
speed,  and  threshold  speed  are  also  input  at  this  time. 

Next,  bins  are  created  for  each  observed  wind  direction  and  peak  speed.  For 
example,  wind  directions  between  0°  and  9°  are  assigned  to  direction  bin  1,  directions 
between  10°  and  19°  are  assigned  to  direction  bin  2,  et  cetera,  until  a  total  of  36 
direction  bins  is  reached.  Wind  speeds  are  handled  in  a  similar  fashion.  Peak  wind 
speeds  between  0  and  4.9  knots  are  assigned  to  speed  bin  1,  speeds  between  5  and 
9.9  knots  are  assigned  to  speed  bin  2,  and  so  on,  until  a  total  of  14  speed  bins  is 
reached. 

Third,  the  number  of  times  each  possible  conditioning  event  (combination  of 
wind  direction  and  peak  speed)  occurred  together  was  totaled.  For  example,  the 
program  counted  the  number  of  times  direction  bin  1  occurred  with  speed  bin  1,  the 
number  of  times  direction  bin  1  occurred  with  speed  bin  2,  et  cetera.  The  resulting 
36  by  14  “conditioning  event  array”  was  saved  for  later  computations.  An  example 
of  a  3  by  3  conditioning  event  array  is  shown  in  Table  10. 

Fourth,  the  program  determines  the  maximum  peak  wind  speed  (the  subse¬ 
quent  event)  which  occurs  for  each  cell  in  the  conditioning  event  array  for  each  of 
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Bins 

1  (0-4.9  kt) 

2  (5-9.9  kt) 

3  (10-14.9  kt) 

1  (0-119°) 

75 

100 

200 

2  (120-239°) 

75 

100 

3  (240-359°) 

50 

100 

100 

Table  10  Sample  3  by  3  conditioning  event  array  for  n  =  1000  observations. 

the  next  eight  hour-long  forecast  periods.  Once  the  maximum  speed  is  found,  it  is 
matched  to  its  corresponding  direction.  Table  11  shows  a  sample  of  the  data  for  one 
hour  and  ten  minutes.  Beginning  at  time  zero  (line  1),  the  direction  is  270°  and  the 
peak  speed  is  8.0  knots.  To  find  the  maximum  peak  speed  over  the  next  hour,  the 
program  examines  the  next  12  observations  (lines  2  through  13).  In  this  case,  the 
maximum  speed  of  13.9  knots  occurs  at  time  4500  (line  10).  Both  the  direction  and 
speed  for  line  10  are  converted  to  their  respective  bins,  and  those  values  are  saved  in 
a  subsequent  event  array.  This  process  is  repeated  for  each  observation.  The  result 
is  a  set  of  504  subsequent  event  arrays,  one  for  each  cell  of  the  conditioning  event 
array.  An  example  of  a  3  by  3  subsequent  event  array  for  the  first  hour  is  shown  in 
Table  12.  This  array  corresponds  to  direction  bin  3  and  speed  bin  2  shown  in  Table 
10. 

Lastly,  each  total  in  the  subsequent  event  array  is  converted  to  a  conditional 
probability  by  dividing  by  the  total  in  the  corresponding  cell  of  the  conditioning  event 
array  (see  Table  13).  These  conditional  probabilities  are  summed  for  each  speed  bin 
greater  than  or  equal  to  the  threshold  bin,  according  to  direction.  For  example,  if 
the  user  was  only  interested  in  maximum  wind  speeds  meeting  or  exceeding  5  knots, 
conditional  probabilities  in  column  1  would  be  discarded,  conditional  probabilities 
in  columns  2  and  3  are  added  together,  and  the  result  is  shown  in  the  last  column 
in  Table  14.  This  step  allows  the  user  to  see  the  total  conditional  probability  of 
meeting  or  exceeding  the  threshold  speed  by  direction,  not  just  the  probability  of 
observing  winds  in  individual  speed  bins.  To  make  the  output  easier  to  interpret 
and  verify,  these  probabilities  are  then  summed  over  all  directions.  In  this  case,  the 
total  conditional  probability  of  meeting  or  exceeding  wind  speeds  of  5  knots  from 
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Line  Number 

Date 

Direction 

Peak  Speed 

1 

99060 

270 

8.0 

2 

99060 

272 

4.9 

3 

272 

4.9 

4 

WHtlfM 

275 

8.0 

5 

99060 

2000 

274 

8.0 

6 

eh 

274 

7 

i— 

273 

9.9 

8 

3500 

273 

9 

EH 

4000 

272 

99060 

4500 

273 

13.9 

11 

99060 

5000 

270 

12.0 

12 

99060 

5500 

269 

12.0 

13 

llHil 

268 

12.0 

14 

268 

16.0 

Table  11  Example  of  data  used  in  determining  maximum  wind  speed  forecasts. 
Direction  is  in  degrees  and  speed  is  in  knots. 


Bins 

1  (0-4.9  kt) 

2  (5-9.9  kt) 

3  (10-14.9  kt) 

1  (0-119°) 

10 

15 

25 

2  (120-239°) 

5 

5 

3  (240-359°) 

5 

5 

Table  12  Sample  3  by  3  subsequent  event  array  corresponding  to  the  100  observed 
cases  shown  in  direction  bin  3  and  speed  bin  2  of  Table  10. 
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Bins 

1  (0-4.9  kt) 

2  (5-9.9  kt) 

3  (10-14.9  kt) 

1  (0-119°) 

351  _  ip 
loo _ ify _ 

—  ”  15 

100  ~~  ,i0 

—  —  25 

100  ”  *ZD 

2  (120-239°) 

EjSifll 

3  (240-359°) 

KHB9I 

E99 

35.  =  ip 
_idq _ ilr _ 

Table  13  The  cells  of  the  subsequent  event  array  are  divided  by  the  corresponding 
cell  of  the  conditional  event  array  to  obtain  conditional  probabilities. 


Bins 

2(5-9. 9  kt) 

3(10-14.9  kt) 

Directional  Probability 

1  (0-119°) 

.15 

.25 

.40 

2  (120-239°) 

.05 

.20 

.25 

3  (240-359°) 

.05 

.10 

.15 

Total  Probability 

.80  or  80% 

Table  14  After  discarding  conditional  probabilities  for  low  wind  speeds,  the  re¬ 
maining  probabilities  are  summed  by  direction. 

any  direction  is  80%.  The  entire  process  described  in  this  section  is  repeated  for 
each  hour  of  the  eight-hour  forecast  period. 

Table  15  contains  a  sample  probability  forecast  for  Tower  1101  at  the  162- 
foot  level  using  data  from  February  1995,  1996,  1997,  and  1998.  Input  conditions 
were:  current  wind  direction  of  225°,  current  peak  speed  of  5  knots,  and  a  threshold 
speed  of  15  knots.  The  output  tells  the  user  the  percent  probability  of  meeting  or 
exceeding  the  given  threshold  speed  from  each  10°  increment  for  each  hour  over  the 
next  eight  horns.  For  example,  for  wind  directions  between  340°  and  349°,  there 
is  a  2.8%  chance  that  winds  will  meet  or  exceed  the  15-knot  threshold  during  the 
fifth,  sixth,  seventh,  and  eighth  hours.  Wind  speeds  can  be  expected  to  meet  or 
exceed  the  threshold  value  from  any  direction  when  the  total  of  any  column  meets 
or  exceeds  30%.  In  this  case,  winds  equal  to  or  in  excess  of  15  knots  are  expected 
to  occur  sometime  during  every  hour  beginning  with  the  fifth  hour. 
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Direction  (degrees) 
(M) 


Direction  (degrees) 

Hr  1 

Hr  2 

Hr  3 

Hr  4 

Hr  5 

Hr  6 

Hr  7 

Hr  8 

0-9 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

10-19 

.0 

.0 

m 

HHI 

.0 

.0 

.0 

.0 

im 

.0 

.0 

30-39 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

50-59 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

60-69 

BHI 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

70-79 

.0 

.0 

.0 

.0 

.0 

.0 

80-89 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

90-99 

.0 

.0 

.0 

.0 

.0 

.0 

100-109 

.0 

.0 

_ 

.0 

.0 

110-119 

.0 

.0 

.0 

.0 

.0 

120-129 

.0 

.0 

.0 

■ 

.0 

130-139 

.0 

.0 

.0 

.7 

2.8 

.0 

.0 

.7 

140-149 

.0 

.0 

.0 

3.4 

4.8 

2.1 

150-159 

.0 

.0 

.0 

.0 

.7 

.0 

160-169 

.0 

.0 

.0 

.0 

.0 

.0 

170-179 

1 

.0 

.0 

.0 

.0 

.7 

6.9 

180-189 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

190-199 

.0 

.0 

.0 

.7 

2.1 

2.8 

.0 

.0 

200-209 

.0 

1.4 

.0 

.0 

2.1 

.7 

.0 

.0 

210-219 

.7 

.7 

1.4 

2.8 

8.3 

6.9 

2.8 

2.8 

220-229 

.0 

2.1 

2.8 

2.1 

6.2 

9.7 

230-239 

USB 

2.1 

3.4 

7.6 

13.8 

8.3 

10.3 

240-249 

.7 

5.5 

2.1 

4.1 

.7 

.7 

4.8 

9.0 

250-259 

.0 

1.4 

4.1 

1.4 

4.1 

3.4 

9.0 

7.6 

260-269 

.0 

2.1 

2.8 

1.4 

4.8 

4.1 

4.1 

270-279 

.0 

.0 

.8 

.0 

.0 

.0 

.0 

280-289 

.0 

.0 

.7 

.0 

.0 

.7 

.7 

.0 

290-299 

1 

.7 

.7 

.0 

.0 

300-309 

191 

.0 

.7 

.0 

.0 

310-319 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

320-329 

.7 

.0 

.0 

.0 

.0 

.0 

.0 

330-339 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

340-349 

.0 

.0 

.0 

2.8 

2.8 

2.8 

2.8 

350-359 

1 

.0 

.0 

2.1 

.0 

2.1 

2.1 

2.1 

Total 

2.1 

11.0 

11.7 

20.7 

35.2 

45.5 

46.9 

57.9 

Table  15  Sample  forecast  for  Tower  1101  at  the  162-foot  level  for  February.  This 
forecast  is  based  on  a  current  direction  of  225  degrees,  current  peak  speed 
of  5  knots,  and  threshold  speed  of  15  knots. 


4.3  Testing  the  Algorithm 

To  ensure  the  program  was  correctly  calculating  the  probability  forecasts,  a 
test  forecast  was  created  using  288  observations  (1  day).  The  forecast  was  then  re¬ 
calculated  by  hand,  and  it  was  determined  that  the  algorithm  was  working  correctly. 

The  algorithm  was  tested  for  accuracy  by  calculating  36  different  conditional 
probability  forecasts  for  each  month  (one  for  each  direction  bin)  using  the  data  from 
the  first  four  winters.  The  fifth  winter  was  used  to  verify  the  forecast.  If  the 
total  conditional  probability  was  equal  to  or  greater  than  30%,  winds  equal  to  or 
greater  than  the  threshold  value  were  expected  to  occur.  If  the  total  conditional 
probability  was  less  than  30%,  strong  winds  were  not  expected.  Contingency  tables 
were  constructed  based  on  these  total  probabilities  and  the  verification  data,  and 
measures  of  accuracy  were  calculated. 

4-4  Results 

Figures  10  through  16  display  measures  of  accuracy,  bias,  and  skill  scores  for 
tower  1101  at  the  162-foot  level.  February  shows  the  best  results  by  far,  with  hit 
rates  above  70%  in  all  time  periods  and  increasing  skill  over  eight  hours.  However, 
the  high  hit  rate  is  due  to  the  relatively  high  number  of  correct  “no”  forecasts  rather 
than  the  number  of  correct  “yes”  forecasts  (see  contingency  table  in  Table  16).  False 
alarm  rates  could  not  be  calculated  for  the  first  three  horns  because  the  “yes”  rows 
in  the  contingency  tables  were  zero.  After  the  first  three  horns,  false  alarm  rates 
range  from  32%  to  51%.  Throughout  the  eight  hours,  critical  success  indices  range 
from  0%  to  24%,  and  probabilities  of  detection  range  from  0%  to  32%.  Bias  ratios 
show  the  algorithm  is  consistently  underforecasting  the  occurrence  of  strong  winds, 
although  they  do  increase  with  time.  Even  though  February’s  scores  axe  the  best 
out  of  all  five  months,  they  still  leave  a  lot  to  be  desired. 

Results  from  November,  December,  January,  and  March  also  show  high  hit 
rates.  However,  false  alarm  rates  are  dramatically  higher,  ranging  from  51%  to 
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Figure  10  Hit  Rate  (H),  Critical  Success  Index  (CSI),  Probability  of  Detection 
(POD),  and  False  Alarm  Rate  (FAR)  for  November,  Tower  1101  at  the 
162-foot  level. 

100%  because  of  the  extremely  high  number  of  incorrect  “yes”  forecasts.  Critical 
success  indices  range  from  0%  to  16%  and  probabilities  of  detection  range  from  0% 
to  27%.  Bias  ratios  reveal  that  the  algorithm  is  underforecasting  strong  winds  in  the 
first  three  hours,  then  begin  to  increase  to  near  one  by  the  eight-hour  point.  The 
only  exception  is  November,  in  which  bias  ratios  show  the  algorithm  overforecasting 
strong  winds.  Heidke  skill  scores  are  negative  or  very  close  to  zero  for  all  four 
months.  Contingency  tables  for  November,  December,  January,  and  March  can  be 
found  in  Appendix  B. 

Where  p-values  could  be  computed  (where  “yes”  forecast  rows  are  not  equal 
to  zero),  p-values  were  found  to  be  zero.  This  shows  that  there  is  an  association 
between  the  observations  and  forecasts,  i.e.,  the  observations  are  dependent  on  the 
forecasts  for  the  0.05  significance  level. 
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Figure  11  Hit  Rate  (H),  Critical  Success  Index  (CSI),  Probability  of  Detection 
(POD),  and  False  Alarm  Rate  (FAR)  for  December,  Tower  1101  at  the 
162-foot  level. 


Figure  12  Hit  Rate  (H),  Critical  Success  Index  (CSI),  Probability  of  Detection 
(POD),  and  False  Alarm  Rate  (FAR)  for  January,  Tower  1101  at  the 
162-foot  level. 
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Hour 


Figure  13  Hit  Rate  (H),  Critical  Success  Index  (CSI),  Probability  of  Detection 
(POD),  and  False  Alarm  Rate  (FAR)  for  February,  Tower  1101  at  the 
162-foot  level. 


Figure  14  Hit  Rate  (H),  Critical  Success  Index  (CSI),  Probability  of  Detection 
(POD),  and  False  Alarm  Rate  (FAR)  for  March,  Tower  1101  at  the 
162-foot  level. 


Observed 


Total 


Hour  1 

Yes 

No 

Forecast 

Yes 

o 

0 

No 

30067 

260201 

P-values 


Div/0 


Hour  2 


Yes 

No 

Yes 

0 

0 

No 

63542 

226726 

Div/0 


Hour  3 


Yes 

No 

Yes 

o 

0 

No 

69749 

220519 

Div/0 


Hour  4 


Yes 

No 

Yes 

5473 

2590 

No 

67230 

214975 

Hour  5 


Yes 

No 

Yes 

12794 

11395 

No 

72298 

193781 

Hour  6 


Yes 

No 

Yes 

12838 

11351 

No 

60847 

205232 

Hour  7 


Yes 

No 

Yes 

10300 

5826 

No 

63350 

210792 

Yes 

No 

Yes 

23741 

24637 

No 

50102 

191788 

Hour  8 


Table  16  Contingency  tables  and  p-values  for  February,  Tower  1101  at  the  162- 
foot  level.  "Total"  column  is  the  total  number  of  observations  for  each 
hour. 
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V.  Conclusion  and  Recommendations 


5. 1  Conclusion 

Conditional  climatology  can  be  a  useful  tool  for  predicting  the  occurrence  of 
strong  winds  at  the  Kennedy  Space  Center  and  Cape  Canaveral  Air  Station.  How¬ 
ever,  as  shown  in  Chapter  4,  accuracy  and  skill  are  severely  limited  in  this  study. 
This  method  of  forecasting  maximum  wind  speeds  is  not  recommended  for  opera¬ 
tional  use. 

5.2  Recommendations  for  Future  Research 

The  results  of  this  study  would  likely  improve  if  more  observations  were  col¬ 
lected  and  included  in  the  dataset.  This  method  could  be  repeated  after  data  from 
at  least  five  more  winters  are  obtained.  Additionally,  different  total  probabilities, 
such  as  10%,  20%,  and  40%  should  be  tested  for  accuracy  and  skill  in  forecasting 
strong  winds.  It  is  possible  that  a  total  probability  other  than  30%  might  be  a  better 
predictor  of  strong  winds.  It  is  also  possible  that  the  “optimal”  total  probability 
could  vary  by  month  or  by  hour. 

A  better  option  for  predicting  strong  winds  at  KSC/CCAS  would  be  to  perform 
a  regression  analysis.  In  addition  to  current  wind  speed  and  direction,  variables 
such  as  temperature  or  pressure  change  with  time,  temperature  or  pressure  gradient, 
wind  speed  or  direction  change  with  time,  gust  spread,  or  wind  speeds  at  upstream 
stations  can  be  included  in  the  study.  Again,  the  accuracy  of  such  a  study  would 
likely  benefit  from  a  larger  dataset,  as  stratifying  the  data  any  further  would  limit 
the  number  of  observed  cases. 
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Appendix  A.  Observed  Frequency  Tables 


Listed  on  the  following  pages  are  observed  frequency  tables  for  the  entire  winter  for 
each  WINDS  tower  and  level  used  in  this  study. 
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0.52 
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0 

Em 
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0.51 

2.32 

1.05 

im 

10.3 

Total 

5.47 

51.2 

22.8 

12.6 

5.94 

1.56 

0 

100 

Table  17  Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
2,  northwest  side,  at  the  90-foot  level,  based  on  40505  observations.  Wind 
speeds  (first  row)  are  in  knots  and  wind  directions  (first  column)  are  in 
degrees. 
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Table  18  Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
2,  southeast  side,  at  the  90-foot  level,  based  on  40513  observations.  Wind 
speeds  (first  row)  are  in  knots  and  wind  directions  (first  column)  are  in 
degrees. 
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0.36  4.70  3.54  2.83 

330-359  I  0.31  3.14  1.47  1.72 


Total  5.1  41.1  25.1  17.1 


Total 
Table  19 
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7.36 

7.79 


5.39 


6.75 

8.76 


12.7 
9.36 
0.14  I  100 


Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
36,  at  the  90-foot  level,  based  on  41772  observations.  Wind  speeds  (first 
row)  are  in  knots  and  wind  directions  (first  column)  are  in  degrees. 
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6.18 
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21.7 
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10.2 

4.67 

1.37 

0.38 

0.37 

100 

Table  20  Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
393  at  the  60-foot  level,  based  on  34062  observations.  Wind  speeds  (first 
row)  are  in  knots  and  wind  directions  (first  column)  are  in  degrees. 
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Table  21  Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
394  at  the  60-foot  level,  based  on  42118  observations.  Wind  speeds  (first 
row)  are  in  knots  and  wind  directions  (first  column)  are  in  degrees. 
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Table  22  Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
397  at  the  60-foot  level,  based  on  42116  observations.  Wind  speeds  (first 
row)  are  in  knots  and  wind  directions  (first  column)  are  in  degrees. 
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Table  23  Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
398  at  the  60-foot  level,  based  on  42197  observations.  Wind  speeds  (first 
row)  are  in  knots  and  wind  directions  (first  column)  are  in  degrees. 
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Table  24  Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
1101  at  the  54-foot  level,  based  on  35425  observations.  Wind  speeds 
(first  row)  are  in  knots  and  wind  directions  (first  column)  are  in  degrees. 
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Table  25  Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
1102  at  the  54-foot  level,  based  on  42105  observations.  Wind  speeds 
(first  row)  are  in  knots  and  wind  directions  (first  column)  are  in  degrees. 
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Table  26  Percent  observed  wind  speed  and  directions  for  the  entire  winter,  Tower 
1102  at  the  162-foot  level,  based  on  42101  observations.  Wind  speeds 
(first  row)  are  in  knots  and  wind  directions  (first  column)  are  in  degrees. 
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Appendix  B.  Contingency  Tables 

Contingency  tables  and  p- values  for  Tower  1101  at  the  162-foot  level  for  November, 
December,  January,  and  March. 
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Observed 


Total 


Hour  1 


Yes 

No 

Yes 

0 

0 

No 

12534 

298506 

P-values 


Div/0 


Hour  2 


Yes 

No 

Yes 

0 

0 

No 

19528 

291512 

Div/0 


Hour  3 


Yes 

No 

Yes 

0 

8640 

No 

27234 

275166 

Hour  4 


Yes 

No 
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0 

34560 

No 

27599 

248881 

Hour  5 


Yes 

No 

Yes 
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60324 

No 

39456 

211104 

Hour  6 


Yes 

No 

Yes 

13 

60467 

No 

39852 

210708 

Hour  7 


Yes 

No 

Yes 

36 

69084 

No 

34357 

207563 

Yes 

No 

Yes 

89 

51751 

No 

33328 

225872 

Hour  8 


Table  27  Contingency  tables  and  p-values  for  November,  Tower  1101  at  the  162- 
foot  level.  "Total"  column  is  the  total  number  of  observations  for  each 
horn. 
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Hour  1 


Forecast 


Observed 


Yes 

No 

Yes 
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Hour  2 
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No 
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0 

0 
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0 

0 
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Div/0 


Div/0 


Div/0 


Div/0 
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Yes 

No 

Yes 
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No 
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321408 


0 


Table  28  Contingency  tables  and  p-values  for  December,  Tower  1101  at  the  162- 
foot  level.  "Total"  column  is  the  total  number  of  observations  for  each 
horn. 
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Hour  3 


Yes 

No 

Yes 

166 

8762 

No 

46767 

265713 

Hour  4 


Yes 

No 

Yes 

166 

8762 

No 

49512 

262968 

Hour  5 


Yes 

No 

Yes 

789 

25995 

No 

35601 

259023 

Hour  6 


Yes 

No 

Yes 

8511 

36129 

No 

48358 

228410 

Hour  7 


Yes 

No 

Yes 

9494 

44074 

No 

60721 

207119 

Yes 

No 

Yes 

20278 

51146 

No 

53947 

196037 

Hour  8 


Table  29  Contingency  tables  and  p-values  for  January,  Tower  1101  at  the  162-foot 
level.  "Total"  column  is  the  total  number  of  observations  for  each  hour. 
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Observed 

Hour  1 

Yes 

No 

Forecast 

Yes 

0 

0 

No 

41762 

279646 

Total 


P-values 


Div/0 


Hour  2 


Yes 

No 

Yes 

0 

0 

No 

64675 

256733 

Div/0 


Hour  3 


Yes 

No 

Yes 

0 

0 

No 

58150 

263258 

Div/0 


Hour  4 


Yes 

No 

Yes 

25 

8903 

No 

52370 

260110 

Hour  5 


Yes 

No 

Yes 

48 

17808 

No 

47554 

255998 

Hour  6 


Yes 

No 

Yes 

14 

17842 

No 

39551 

264001 

Hour  7 


Yes 

No 

Yes 

31 

17825 

No 

49658 

253894 

Yes 

No 

Yes 

8969 

35671 

No 

43553 

233215 

Hour  8 


Table  30  Contingency  tables  and  p-values  for  March,  Tower  1101  at  the  162-foot 
level.  "Total"  column  is  the  total  number  of  observations  for  each  hour. 
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